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Abstract 



The first version of a computer program eett6f for calculating cross sections of e — > 
6 fermions processes relevant for a tt-pair production and decay at centre of mass energies 
typical for linear colliders is presented. eett6f v . 1.0 allows for calculating both the total 
and differential cross sections at tree level of the Standard Model (SM). The program can be 
used as the Monte Carlo generator of unweighted events as well. 
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PROGRAM SUMMARY 



Title of program: eett6f 
Version: 1.0 (April 2002) 
Catalogue identifier: FJA0815b 

Program obtainable from: CPC Program Library or on request by e-mail from the author 
Licensing provisions: none 
Computers: all 

Operating systems: Unix/Linux 
Programming language used: FORTRAN 90 
Memory required to execute with typical data: 1 Mb 
No. of bits in a word: 32 

No. of bytes in distributed program, including test data, etc.: 400 Kb 
Distribution format: ASCII 

Keywords: e + e~ annihilation, SM, lowest order six-fermion reactions, tt production and decay 
Nature of physical problem 

Description of e + e~ — > 6 fermions processes relevant for a tt-pair production and decay at 
centre of mass energies typical for linear colliders; all six fermion reactions containing a b and 
b quarks and four other fermions of different flavours with a complete set of the Feynman 
diagrams in the lowest order of SM. 

Method of solution 

Matrix elements are calculated numerically with the helicity amplitude method. Constant 
widths of unstable particles are implemented by modifying mass parameters in corresponding 
propagators. The phase space integration is performed numerically utilizing a multi-channel 
Monte Carlo method. 

Restrictions on complexity of the problem 

Reactions containing fermions of the same flavour are not treated. No higher order effects arc 
taken into account, except for assuming the fine structure constant and the strong coupling 
at appropriate scale and partial summation of the one particle irreducible loop corrections by 
introducing fixed widths of unstable particles. 

Typical running time 

The running time depends strongly on a desired precision of the result. The results of the 
appended test run have been obtained on a 800 MHz Pentium III processor with the use of 
Absoft FORTRAN 90 compiler in 490 seconds. In order to obtain a precision level below one per 
mille a few million calls to the integrand are required. This results in a typical running time 
of several hours. The running time becomes much shorter for approximated cross sections. 



LONG WRITE-UP 



1 Introduction 

Precise measurements of the top quark properties and interactions are planned at TESLA 
UK and will most certainly belong to the research program of any future e + e~ collider [Q. 
The measurements should be confronted with theoretical predictions matching the same high 
precision level of a few per mille. It is obvious that in order to reach that high precision 
it is mandatory to include radiative corrections and, as the measurements of some top quark 
physical properties will be carried out at high energies, much above the ti threshold, it is crucial 
to know off-resonance background contributions to any specific 6 fermion decay channel and 
to estimate the effects related to the off-shellness of the tt-pair. 

In this article, a technical documentation of a numerical program, eett6f , is presented which 
allows for computer simulation of the 6 fermion reactions which are relevant for the top quark 
pair production and its decay into a specific 6 fermion final state 

e + ( Pl )+e-(p 2 ) - b{p 3 ) + f 1 {p,) + f[{p 5 )+b(p 6 )+f 2 (p 7 )+h(p 8 ), (1) 

where fi = z/^, z/ r , u, c, f 2 = [i~ ,r~ , d, s, and /(, f 2 are the corresponding weak isospin partners 
of fi, f 2 , f[ = fi~ ,r~ , d, s, f 2 = z/^, z/ r , u, c, and the particle four momenta have been indicated 
in the parentheses. In the present version, v . 1 . 0, of the program, reactions ([!]) are treated in 
the lowest order of SM. Moreover, for the sake of simplicity, it is assumed that the actual values 
of fi and f'2 in ([[]) are different from each other, and that neither f[ nor f 2 is an electron. The 
physics contents of the program is described in B, where the details on the physical model, a 
method of calculation and physical results obtained with eett6f vl.O have been discussed. 

A similar analysis of the 6 fermion processes relevant for a ti production in e + e~ annihilation 
have been performed in 0], ||, where semileptonic channels of reaction ([!]) have been studied, 
and in 0, where purely hadronic channels of (|1]) have been analysed. Moreover, irreducible 
QCD background to top searches in semileptonic channels of ([!]) has been discussed in 0. After 
this work has already been completed an exstensive study of six fermion reactions, including (Q) 
among others, in the massless fermion limit has appeared ||. Using the same input parameters 
and separation cuts as in ||, the present program gives 5.8160(32) fb and 17.223(15) fb for 
the cross sections of e + e" — > bv^fi + bT~v T and e + e~ — > budb\i~v^ respectively, at CMS energy 
of 500 GeV. These results are 2-3 standard deviations off the corresponding results of ||. 
This difference can most probably be traced back to somewhat different implementation of 
the top quark width in the top quark propagator. In the present work, the top quark width 
is introduced both in the denominator and in the numerator of the top propagator while, 
in [|], it seems to be introduced in the denominator only. The invariant mass and angular 
distributions of the bud quark triple of e + e~ — > budb^v^ at s 1//2 = 500 GeV shown in M are 
nicely reproduced by eett6f within accuracy of the plots. 

The main advantage of the present program is that it allows for taking into account both the 
electroweak and QCD lowest order contributions. As light fermion masses are not neglected, 
the cross sections can be calculated without any kinematical cuts. Moreover, a number of 
options have been implemented in the program which make possible calculation of the cross 



sections while switching on and off different subsets of the Feynman diagrams. It is possible 
to calculate cross sections in two different narrow width approximations discussed below, too. 



Besides of a possibility of taking into account solely the electroweak contributions, or switching 
off the Higgs boson exchange diagrams, the program allows also for a simplified treatment of 
reaction (|1|) by utilizing a few different approximations: the double resonance approximation 
for W bosons 

e + bi) + e-(p 2 ) - b(p 3 ) + W + *( P 4 5 ) + b( Pfi ) + W-*(p 78 ) (2) 

- b(Ps) + fM) + f{{Pn) + KP6)+MP7)+h(ps), (3) 

where only those 61 Feynman diagrams are taken into account which contribute to e + e~ — > 
bW + bW~ and the W bosons are considered as being off-mass- shell and the double resonance 
approximation for a t- and t-quark 

e + {pi) + e-{p 2 ) -> t*{p U5 )+i*{p 678 ) (4) 

- b{p 3 ) + f 1 (p 4 ) + f{{p 5 ) + b(p 6 ) + f 2 (p 7 )+n(p 8 ), (5) 

with only two 'signal' diagrams contributing. The intermediate state momenta are P45 = 
P4 + Ps, P78 = Pi + Ps in Eq. (|) and p 345 = p 3 + p 4 + p B , p e78 = p 6 + p 7 + p 8 in Eq. (@). 
The narrow width approximation for the W bosons (top quarks) is obtained by replacing the 
intermediate state W bosons (top quarks) in Eq. (@) (in Eq. (§])) with on mass shell W's 
(top quarks) and by multiplying the cross section of the on mass shell reaction (g) ((§)) with 
the branching ratios corresponding to final state of reaction (p]). The necessary partial decay 
width of the W and t are calculated according to the lowest order SM, and the fermion masses 
are neglected in the three body top decay with. 

The necessary matrix elements of reactions ([!]), (fj) and (Q) are calculated with the helicity 



amplitude method described in and [10) and phase space integrations are performed with 
the Monte Carlo (MC) method. More details on the multi-channel MC algorithm (see e.g. 
|TT|), on which the MC integration and event generation is based, are given in the next section. 



2 Monte Carlo integration and event generation 

In order to improve the convergence of the MC integration the most relevant peaks of the 
matrix element squared related to the Breit-Wigner shape of the W, Z, Higgs and top quark 
resonances as well as to the exchange of a massless photon or gluon have to be mapped away. 
As it is not possible to find out a single parametrization of the n- dimensional phase space 
which would allow to cover the whole resonance structure of the integrand, it is necessary to 
utilize a multi-channel MC approach. In the multi-channel MC approach, the random numbers 
< Xj < 1, j = 1, ...,n, and hence four momenta of the final state particles, are generated 
according to one of iV different probability densities fi(x), i — 1, N, x = (x\, ...,x n ). The 
distribution fi(x), which accounts for several different peaks of the integrand, is selected with 
probability a^. All the distributions fi(x) are then combined in a single probability distribution 
f(x) which should cover possibly all the peaks of the integrand, or at least the most relevant 
ones 

N N 

f( x ) = J2 a ifc( x )i with J2 a * = 1 - ( 6 ) 

i=i i=i 



The second condition of Eq. (|5|) guaranties that, if every distribution fi(x) is normalized to 
unity, then the combined distribution f(x) is normalized to unity, too. Obviously, all the 
distributions fi(x) and weights a iy i — 1, N, are to be non negative. 

The basic phase space parametrizations which are used in the program are listed below. The 
6 particle phase space of reaction (|1|) is parametrized in 3 different ways: 

d u Lips = (2ir)~ 14 dPS 2 (s, s 345 , s 678 ) dPS 2 (s 345 , ml, s 45 ) dPS 2 (s 678 , m\, s 78 j 

x dPS 2 (s 45 , m\, mfj dPS 2 (s 78 , m 7 , mfj ds 345 ds 678 ds 45 ds 78 , (7) 

d u Lips = (2iry u dPS 2 (s, s 34 , s 5678 ) dPS 2 (s 5678 , s 56 , s 78 ) dPS 2 (s 34 , m\, m\ 

x dP^ (s 5 e, ml, mfj dPS 2 (s 78 , m 7 , mfj ds 34 ds 5678 ds 56 ds 78 (8) 

and 

d u Lips = (27r)~ 14 dPS 2 (s, m\, s 45678 ) dPS'2 45 678, s 45, S67s) dPS 2 (s 678 , m\, s 7 fj 

x dPS 2 [S45, ml, mfj dPS 2 (s 78 , m 7 , mfj ds 45678 ds 45 ds 678 ds 78 . (9) 

The 4 particle phase space of reaction (fj) in the narrow W width approximation is parametrized 
in 2 different ways: 



d Lips = (2tt) dPS 2 (s, m 3 , s 456 J dPS 2 [s 456 , m 4 , s 56 J dPS 2 [s m , m 5 , mfj ds 456 ds 56 (10) 
and 

d 8 Lips = (27r) _8 dPS 2 (s, s 34 , s 56 ) dPS 2 (s 34 , m\, m\, ) dPS 2 (s 56 , m\, mfj ds 34 ds 56 . (11) 

In Eqs. (|7|-4TlD, s^k... = (Pi+Pj +Pk + ■■■) 2 , hj, k = 3, 8, and dPS 2 (s, s', s") is a two particle 
(subsystem) phase space element defined by 

dPS 2 (s, s', s") =5 A (p-p'- p") f^fPl = Mld^', (12) 
v ' ' ' KF F F ' 2E' 2E" A^/7s ' v ; 

where p ' is the momentum and Q' is the solid angle of one of the particles (subsystems) in 
the relative centre of mass system, p ' + p " = 0. Making use of the rotational symmetry with 
respect to the e + e~ beam axis reduces the dimension of the phase space elements to 13 in 
Eqs. ®-(D and to 7 in Eqs. (0) and (0). 



Parametrizations ([7HTTD are used with different permutations of external particle momenta 
such that invariants possibly correspond to the virtuality of propagators of the gauge 
bosons, Higgs boson and/or top quarks. The invariants are then transformed to the interval 
[0, 1] by performing mappings smoothing out peaks related to the propagators. All other 
integration variables of Eqs. (|7HTTD, not related to peaks of the matrix element squared, are 
also transformed to the interval [0, 1] with a simple mapping defined in the following way. 
Let yj G [cj,bj] be a variable not related to a peak in the matrix element squared and Xj be 
random variable uniformly distributed in the interval [0, 1], then the relation 



Vj = i h i ~ '•,) x 3 + 'j 



(13) 



defines the necessary mapping. After having all the integration variables transformed to the 
interval [0,1], any of phase space parametrizations (0^0) for reaction ([[]), or (|10D , ( pTTj ) for 
reaction @ in the narrow width approximation, for a given permutation of particle momenta, 
can be directly associated with the probability density fi(x), which is referred to as a kinemat- 
ical channel. Altogether 59 (21) different channels are used in order to integrate the matrix 
element squared of reaction ([]]) (reaction (Q) in the narrow width approximation). 

The weights a% are calculated in the initial scanning run that starts with all the weights equal 
to each other. The weights aj, % = 1, N, for the actual run are calculated as ratios 

JV 

Cli = (Ti/^CTj, (14) 

i=i 

where Oj denotes the cross section obtained with the i-th kinematical channel in the initial 
scan. 

Calculation of the cross section of reaction (||) in the narrow top width approximation does 
not require the multi-channel MC approach and can be performed with the single phase space 
parametrization ([12]). 



3 Description of the program 

eett6f is a package written in FORTRAN 90. It consists of 50 files including a makefile. They 
are stored in one working directory. The user should specify the physical input parameters in 
module inprms.f and select a number of options in the main program csee6f .f. 

3.1 Program input 

The default values of the input parameters and options used in the program are those specified 
below. 

3.1.1 Physical parameters 

The initial physical parameters to be specified in module inprms . f are the following. 

Gauge boson masses and widths (GeV): 

mw=80.419 GeV, the W mass, 
gamw=2 . 12 GeV, the W width, 
mz=91 . 1882 GeV, the Z mass, 
gamz=2 . 4952 GeV, the Z width. 

The electroweak (EW) mixing parameter sw2 is then calculated from 

sw2 = l-mw2/mz2 
utilizing mw2 = mw**2, mz2 = mz**2 in the fixed width scheme, and 



mw2 = mw**2 - i*mw*gamw, mz2 = mz**2 - i*mz*gamz in the complex mass scheme. 
ralpO = 137.03599976, an inverse of the fine structure constant in the Thomson limit, 
gmu=l . 16639x 10 -5 GeV -2 , the Fermi coupling constant, 
alphas=0. 1185, the strong coupling constant at mz. 

An inverse of the fine structure constant at mw, ralpw, is then calculated from 

ralpw=4 . 44288293815837/ (2*sw2*gmu*mw**2) . 
mh=115.0 GeV, the Higgs boson mass, 
gamh=0 GeV, the Higgs boson width. 

If gamh = 0, then the Higgs width is calculated according to the lowest order of SM. 
Fermion masses and widths: 

me = 0.510998902 MeV, game = MeV, for an electron, 
mmu = 105.658357 MeV, gammu = MeV, for a muon, 
mtau=l .77703 GeV, gamtau = GeV, for a lepton tau, 
mu = 5 MeV, gamu = MeV, for an up quark, 
md=9 MeV, gamd = MeV, for a down quark, 
mc=1.3 GeV, game = GeV, for a charm quark, 
ms=150 MeV, gams = MeV, for a strange quark, 
mt=174.3 GeV, gamt = 1.5 GeV, for a top quark, 
mb=4 . 4 GeV, gamb = GeV, for a bottom quark. 

ncol=3, the number of colours, 

convc=0.389379292xl0 12 fb GeV 2 , a conversion constant. 
3.1.2 Main options 

The following main options should be selected in the main program csee6f . f . 

The number of different centre of mass (CMS) energies ne 

ne = 1. Recommended if the unweighted event are to be generated. 

The actual values of the CMS energies in the array aecm of size ne: 

aecm=(/500.d0/). 

The final state of ([[]) by selecting a decay mode of the intermediate W + , W~ bosons, iwp, 
iwm = 1 or 2, where 1 corresponds to a leptonic, and 2 to a hadronic decay mode, and by 
specifying the family index if p , if m = 1 , 2 , 3 of the fermion pair resulting from the W + or W~ 
decay. For example, the reaction e + e~ — > bv^^bdu corresponds to W + decaying leptonically 
and W~ decaying hadronically and the corresponding flags are: 

iwp = 1 
ifp = 2 
iwm = 2 



Whether or not to calculate the Born cross section utilizing the full 6 particle kinematics, 
ibom = 1 (yes) / else (no), with ncallO calls to the integrand 

iborn = 1 

ncallO = 20000. Recommended No of calls is a few millions. 
Generate the unweighted events or not, imc = 1 (yes) /else (no)? 
imc = . 

No standard event record is used. If imc = 1, then the final state particle momenta of the 
accepted unweighted events are printed in the output. 

Scan the Born cross section with nscanO calls, iscanO = 1 (yes) /else (no), in order to find 
the dominant kinematical channels, adjust integration weights and find out the maximum 
value of the cross section 

iscanO = 1. This option is strongly recommended. 
nscanO = 200. About 1 thousand is recommended. 

Calculate an approximate cross section of (§), iwwa = 1 (yes) / else (no)? 
iwwa = 1 

ncallww = 20000. Recommended value is several hundred thousands. 

Scan the Born cross section in the narrow W- width approximation with nscanww calls, iscanww 
= l(yes)/else(no)? 

iscanww = 1. Recommended. 

nscanww = 200. Recommended value is about 1 thousand. 

Calculate an approximate cross section of (|j), itopa = 1 (yes) /else (no)? 

itopa = 1 
ncalla = 20000. 

Include the Higgs boson exchange, ihiggs = 1 (yes) /else (no)? 
ihiggs = 1 

Calculate the electroweak contributions only, iew = 1 (yes) /else (electroweak + QCD 
Feynman diagrams)? 

iew = 0. 

Calculate the ti resonant diagrams only, itt = 1 (yes) /else (all the diagrams)? 
itt = 0. 

Calculate the double W resonant diagrams only, iww = 1 (yes) /else (all the diagrams)? 
iww = 0. 

Choose the scheme: ischeme = 1 (complex mass scheme) /else (fixed width scheme) 
ischeme = 1. 

If ischeme=l, then should alpha_W be complex (iaplw=l) or real (ialpw=0)? 
ialpw=l . 



3.1.3 Auxiliary options 



Impose cuts, icuts = 1 (yes) /else (no)? 
icuts = 

If icuts = 1, then specify the kinematical cuts, e.g.: 

ctlb = cos(10*pi/180) - cosine of the charged lepton-beam angle, 

ctqb = cos(5*pi/180) - cosine of the quark-beam angle, 

ctll = 1 cosine of the charged lepton-charged lepton angle, 

ctlq = cos(5*pi/180) - cosine of the charged lepton-quark angle, 

ecutl = 1 GeV - minimum charged lepton energy, 

ecutq = 3 GeV - minimum quark energy, 

mqq = 10 GeV - minimum invariant mass of a quark pair, 

mil = GeV - minimum invariant mass of a charged lepton pair, 

Calculate distributions, idist = l(yes)/0(no)? 

idis=0. 

If idis = 1, then specify parameters of the distributions: 
xmin = (/ xll, xl2, xl3, xl4/) - lower bounds, 
xmax = (/ xul, xu2, xu3, xu4/) - upper bounds, 

nbs = (/ nl, n2, n3, n4/) - the corresponding numbers of bins in each distribution. 

Constants xli, xui should be of type real (8) and ni of type integer. The number of de- 
sired distributions and the maximum number of bins, nbmax = max{nl, n2, n3, n4} should 
be specified in module distribs.f. 

The maximum of the fully differential cross section crmax, relevant only if iscanO = or 
iscanww = 0, 

crmax=1000. 



3.2 Routines of eett6f 

The main program csee6f . f , each subroutine, function or module are located in a file named 
the same way as the routine itself, except for 3 functions: srr, src and sec located in a 
single file dotprod.f. The program consists of the following routines. 

• SUBROUTINE boost - returns a four vector boosted to the Lorentz frame of velocity -v. 

• SUBROUTINE carlos - the MC integration routine. 

• SUBROUTINE couplsmb - returns the SM couplings. 



FUNCTION cross - calculates the cross section of (Ij), approximated cross sections of 
and the cross section of (|2|) in the narrow P^-width approximation. 

FUNCTION crosstopa - calculates the cross section of (f|) in the narrow top width ap- 
proximation. 

The MAIN PROGRAM csee6f . 

MODULE distribs - contains parameters of distributions. 

MODULE drivec - contains driving flags and some kinematical variables. 

SUBROUTINE eebwbw - returns the squared matrix element averaged over initial spins 
and summed over final spins and colours of e + e~ — > bW + bW~ . 

SUBROUTINE eee - returns a contraction of a triple gauge boson coupling with three 
complex four vectors. 

SUBROUTINE eeee - returns a contraction of a quartic gauge boson coupling with four 
complex four vectors. 

SUBROUTINE eeff 1 - returns the squared matrix element averaged over initial spins and 
summed over final spins and colours of e + e~ — > ff. 

SUBROUTINE eeh - returns a contraction of the Higgs-gauge boson coupling with two 
complex four vectors. 

SUBROUTINE eett6f - returns the squared matrix element averaged over initial spins 
and summed over final spins and colours of reaction ([D) and approximations 

SUBROUTINES eev, eve and vee - return contractions of a triple gauge boson coupling 
with two polarization vectors leaving, respectively, the third, second and first Lorentz 
index uncontracted. 

SUBROUTINE fhfb - returns a set complex scalars: scal(\i, A 2 ) = gsi2UiU 2 A F , where Ui, 
U2 are fermion spinors, Ap is a scalar boson propagator and gsn is a coupling of the 
scalar boson to fermion 1 and 2. 

SUBROUTINE f vf a - returns a set of four vectors £ M (Ai, A 2 ) = u^X^^g^ P_ + gfr ] P + ) 
XM 2 (A 2 )-Dy\ where u±, w 2 are fermion spinors, P± = (1 ±7 5 )/2, are chirality projectors, 
g^ are the fermion-gauge boson couplings of definite chirality and Dy 1 is the propagator 
of a gauge boson V in arbitrary linear gauge. 

SUBROUTINE f ef - returns a set of matrix elements mat(\i, A, A 2 ) = Ui(\i)$(\)(g v ^P_ + 
gy^P + )w 2 (A 2 ), where notation is the same as in fvf a. 

SUBROUTINE impcuts - imposes cuts on the kinematical variables of process (jl]). 
MODULE inprms - contains initial input parameters. 
MODULE kincuts - contains kinematical cuts. 

SUBROUTINE kineel212 - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 6 process in the CMS. The phase space parametrization used is that of 
Eq. ©. 



SUBROUTINE kinee6f - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 6 process in the CMS. The phase space parametrization used is that of 
Eq. ©• 

SUBROUTINE kineeff 1 - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 2 process in the CMS. The phase space is parametrized according to 
Eq. 

SUBROUTINE kineef f f f - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 4 process in the CMS. The phase space is parametrized according to 
Eq. (0). 

SUBROUTINE kineetbf f - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 4 process in the CMS. The phase space is parametrized according to 
Eq. (13). 

SUBROUTINE kineett6f - returns the four momenta, phase space normalization and flux 
factor for a 2 — > 6 process in the CMS. The phase space parametrization used is that of 
Eq. (§. 

SUBROUTINE kinf f - returns the final state four momenta and phase space normalization 
for a 2 — > 2 process in the CMS. 

FUNCTION lamsq - the kinematic lambda function, X(y/x, y/y, \/z)- 
MODULE mathprms - contains necessary arithmetical constants. 

SUBROUTINE parfix - returns parameters of a specific process ([I]) and the branching 
ratios for the narrow width approximations; initializes weights for the Monte Carlo 
integration. 

MODULE parproc - contains parameters for a specific process. 

SUBROUTINE peu - returns a set of generalized spinors: u(X, Ai) = SF$(X)(g v P_ + 
gy P + )ui(Xi), where Sp is a Feynman propagator of an internal fermion and the re- 
maining notation is the same as in fvf a. 

SUBROUTINE psneel212 - returns a phase space normalization as the one of kineel212 
for a given set of external particle momenta. 

SUBROUTINE psnee6f - returns a phase space normalization as the one of kinee6f for a 
given set of external particle momenta. 

SUBROUTINE psneef f f f - returns a phase space normalization as the one of kineef f f f 
for a given set of external particle momenta. 

SUBROUTINE psneetbf f - returns a phase space normalization as the one of kineetbf f 
for a given set of external particle momenta. 

SUBROUTINE psneett6f - returns a phase space normalization as the one kineett6f for 
a given set of external particle momenta. 

SUBROUTINE recpol - returns real polarization vectors of a vector boson in the rectan- 
gular basis. 



• FUNCTION sec (contained in dotprod.f) - returns the Minkowski dot product of two 
complex four vectors. 

• SUBROUTINE spheric - returns spherical components ps of a four vector p M , ps= 

( | p| , cos 9, sin 9, cos </>, sin <p) , with 9 and being a polar and azimuthal angles of momen- 
tum p. 

• SUBROUTINE spine - returns the contractions: p°I — pa and p°I + p • er, where p 11 = 
(p°, p) is a complex four vector, / is the 2x2 unit matrix and a are the Pauli matrices. 

• SUBROUTINE spinornew - returns helicity spinors in the Weyl representation. 

• SUBROUTINE spinr - returns the same contractions as in spine for a real four vector ft 1 . 

• FUNCTION sre (contained in dotprod.f) - returns the Minkowski dot product of a real 
and a complex four vector. 

• FUNCTION srr (contained in dotprod.f) - returns the Minkowski dot product of two 
real four vectors. 

• SUBROUTINE uep - returns a set of generalized spinors: w(Ai,A) = Ui(\i)$(\)(g^ P_ + 
g^P^Sp, where notation is the same as in peu. 

• SUBROUTINE wwidth - returns the partial width of the M^-boson averaged over initial 
spins and summed over final spins, and colours for hadronic decay modes. 



3.3 Run output 

A sample of the listing of the test run output is given in Appendix A. It contains a specification 
of the process, information on the scheme choice, and values of the relevant physical parameters 
in the very beginning. Then, for each value of the CMS energy and for each cross section 
calculated, results of the initial scan and resulting weights for the actual run are printed, if 
iscan = 1. Sometimes the weights calculated according to Eq. fll4|) do not add up exactly 
to 1. Then, the last weight is changed slightly in order the second condition of Eq. (^) to be 
fulfilled. This change is supposed to be completely irrelevant numerically, but the message 
about it is printed. The final result for the total cross section is called Integral. Its value 
is printed together with the standard deviation and the actual number of calls used in the 
calculation. Finally, information on events acceptance efficiency is given, that means a fraction 
of accepted weight 1 events. If the event generation option is switched on, i.e. imc = 1, then 
the corresponding number of unwighted events will be printed out as collections of final state 
particle four momenta. 

Whenever a maximum value of the cross section initially assigned in csee6f . f or found in 
the result of the initial scan is overflown by more than a factor 1.5, a corresponding message 
informing about it is printed. The integration is still valid, however, if the program is run 
as an event generator, it should be rerun. How to proceed in this situation is described in 
the next section. In the very end, all the calculated total cross sections together with the 
corresponding standard deviations are collected in the tabular form. 



4 Use of the program 



Up to now the program have been run only on Unix/Linux platforms. In order to run the 
program, the user should select a specific name of a FORTRAN 90/95 compiler, desired options 
and the name of the output file in the makefile. The output file is called test at present. 
The program can be then compiled, linked and run by executing a single command 

make test. 

The results of the test run should reproduce those contained in file testrun. 

The program can be run as the MC event generator of unweighted events by selecting 
imc = 1 

in csee6f .f. It is then recommended to run the program for a single CMS energy and 
to perform a scan with a relatively large number of calls, nscanO or nscanww. Attention 
should be paid to possible messages informing about updates of the maximum weight. In this 
case the program should be rerun, however, this time without the initial scan. The initial 
integration weights awO or awa in parfix.f and the maximum weight crmax in csee6f .f 
should be updated according to the results of the prior scan. At present the efficiency of 
events acceptance is relatively low. However, in view of a relatively fast performance of the 
program, this should not be a serious limitation for the user. 



5 Outlook 



The following improvements of the program are envisaged in the near future. First of all, a 
new version is being prepared which will allow for calculations of reactions involving the same 
external fermion flavours. Another step will be an inclusion of the most relevant factorizable 
higher order effects: electroweak and QCD corrections to the top pair production, e + e~ — > ti, 
to the t and i quark decays, as well as to the W + and W~ decays. 
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Appendix 
A Listing of test run output 



Approaching the process: 

e+ e- -> b vm mu+ b d u 

with the program 'eett6f v. 1.0' by 

Karol Kolodziej 
University of Silesia 
Katowice, Poland 
e-mail: kolodzie@us.edu.pl 



Complex mass scheme 

with sw2 = ( 0.222271237327414,-0.000778147496385) 

and complex l/alpha_W = (132.386398893114400, 3.953796222427256) 

Gauge boson masses and widths: W, Gm_W, Z, Gm_Z 
80.41900 2.12000 91.18820 2.49520 

Higgs boson mass and width: H, Gm_H 
0.11500E+03 0.49657E-02 

1/alphaO = 137.0359998 alphas = 0.1185000 Gjiu = 0.116639E-04 

t-quark mass and width: 
174.30 1.50000 



External fermion masses: e- b vm mu+ d u 



0.51099890E-03 . 44000000E+01 . OOOOOOOOE+00 . 10565836E+00 
0.90000000E-02 . 50000000E-02 



E_cm = 500.000 GeV 



All cross sections in femtobarns (fb) 



The Born cross section with full 6 particle kinematics: 



Results of the initial scan: 

Kinematics Weight Integral St. dev. No of calls 

1 0.17241379E-01 . 77943026E+01 . 53027331E+00 200 

2 0.17241379E-01 . 70433588E+00 . 14557680E+00 200 

3 0.17241379E-01 . 46787757E+00 . 14942639E+00 200 

4 0.17241379E-01 . 40042444E-01 . 17572985E-01 200 

5 0.17241379E-01 . 61761884E-01 . 17401289E-01 200 

6 0.17241379E-01 . 70179588E-02 . 50890185E-02 200 

7 0.17241379E-01 . 18918195E-04 . 14788334E-04 200 

8 0.17241379E-01 . 59135660E-02 . 26123591E-02 200 

9 0.17241379E-01 . 65946999E-05 . 37561359E-05 200 



10 


0. 
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■01 


0. 


.71059323E- 


02 


0. 
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200 


11 
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0. 
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01 
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01 
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02 
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01 


0. 
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•01 


0. 
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Born scanned with 11565 calls = . 19039197E+02 

The maximum value of the cross section= . 27872 147E+04 

shifted to 0.41808220E+04 



Weight 58 = . 4148134056762879E-03 rescaled to . 4148134056756181E-03 
Calculating the Born cross section with the weights: 

1 0.40938191E+00 . 40938191E+00 

2 0.36993992E-01 . 44637590E+00 

3 0.24574439E-01 . 47095034E+00 

4 0.21031583E-02 . 47305350E+00 

5 0.32439333E-02 . 47629743E+00 

6 0.36860583E-03 . 47666604E+00 

7 0.99364459E-06 . 47666703E+00 

8 0.31059956E-03 . 47697763E+00 

9 0.34637490E-06 0.47697798E+00 

10 0.37322648E-03 . 47735120E+00 

11 0.64769547E-02 . 48382816E+00 

12 0.26653539E-02 . 48649351E+00 

13 0.27491290E-02 . 48924264E+00 

14 0.87699297E-04 . 48933034E+00 

15 0.11411925E-03 . 48944446E+00 

16 0.52889841E-02 . 49473344E+00 

17 0.29107678E-02 0.49764421E+00 

18 0.25048428E-04 . 49766926E+00 

19 0.64448631E-07 . 49766932E+00 

20 0.30273882E-03 .49797206E+00 

21 0.57506089E-05 . 49797781E+00 



99 






-OA 




AQ7QQPM RF+nn 


9^ 


n 




-04 


n 




94 


n 




uu 


n 




9R 


n 


. OS i i OSOOL 


uo 


n 




9fi 
zu 


n 
w . 


.21284769E- 


•02 


n 


RnDR4-9DRF+nn 

. JU\JUtZ\JOJ_j t UU 


97 


n 


.25393213E- 


•04 


n 


R00Rfi74.7F+nn 


9R 
zo 




.92912076E- 


■07 




R0nRfi7R7F+00 


9Q 


n 


. 13101157E- 


•04 


n 


R00Rsnfi7F+nn 


ou 


n 


.35784246E+00 






O X 


n 


.42088930E- 


•01 




QOOR1 onfiF+nn 

. yuUOizuoijTuu 


•^9 
oz 


n 


, 19449451E- 


•01 




Q1 QQR^ R1 F+OO 




n 


, 14686352E- 


•02 


n 


Q91 A^m AF+nn 




n 


.65033652E- 


•04 


n 


Q91 AQR1 RF+00 




n 


20278362E- 


06 


n 


Q91 AQR^RF+OO 




n 


. 11861960E- 


•04 


n 


Q91 Rn794.F+00 


^7 


n 


.53899508E- 


•07 


n 


Q91 R079QF+00 


oo 




.31457348E- 


05 


n 


Q91 R1 OAzlF+nn 




n 


.62008511E- 


06 


n 


Q91 R1 1 OfiF+00 


zin 


n 


. 19686903E- 


•03 


n 


Q91 707Q^F+00 




n 


.36746070E- 


03 


n 


Q9907R^QF+00 


zl.9 
^z 




. 13822779E- 


•04 


n 


Q990SQ91 F+00 


A/3, 


n 


.78687301E- 


06 


n 


Q99nQnnnF+nn 


Azl 


n 


. 19213421E- 


•03 


n 


Q999S91 ^F+00 


*±o 


n 


.34504971E- 


06 


n 


Q999894PF+00 




n 


.31137396E- 


02 


n 


Q9R^Qfi99F+00 

. JZUOJUZZEiTvU 


47 


n 


. 18365550E- 


•03 


n 


Q9RR7QS7F+00 


*XO 


n 


.20729799E- 


•02 


n 


Q97fiR9RRF+On 


zlQ 


n 


. 17593374E- 


•03 


n 


Q9789a7QF+nn 
. 3z r ozo f yHjTVju 


RO 
ou 




. 34378833E- 


01 


n 


Qfi9907fi9F+00 


O X 


n 


.61849819E- 


03 


n 


Qfi9S9fi1 9F+00 


R9 


n 


.65075477E- 


03 


n 


Qfi^4.7fiR7F+nO 


RQ 


n 


. 10004530E- 


02 


n 


QfiZLzl77^^F+nn 


V.A 

o^± 


n 


. 15495510E- 


•02 


n 


Qfifin9fiR«F+nn 

. jODUZDOOLtWU 


RR 

oo 


n 


32832332E- 


•01 


n 


QQSSRQ91 F+00 

. 3jOOUi7Z1Ij t uu 


56 


. 


43282221E- 


•03 


. 


. 99929203E+00 


57 


0. 


29315502E- 


•03 


0. 


.99958519E+00 


58 


0. 


.41481341E- 


•03 


0. 


, 10000000E+01 


59 


0. 


.00000000E+00 


0. 


.00000000E+00 


60 


0. 


.00000000E+00 


0. 


.00000000E+00 



Final result : 






Integral = 0.201003E+02 Standard dev. 


= 0.27E+00 


with 19980 calls 


Efficiency of the events acceptance = 


= 0.410E-02 




The narrow W-width approximation, e+ 


e > b W+ b 


W- -> 6f: 



Results of the initial scan: 

Kinematics Weight Integral St. dev. No of calls 



1 0.47619048E-01 . 23618560E+00 . 10375292E+00 200 

2 0.47619048E-01 . 25574899E+00 . 11393103E+00 200 

3 0.47619048E-01 . 90318927E-02 . 27624478E-02 200 

4 0.47619048E-01 . 14178256E+00 . 63376575E-01 200 

5 0.47619048E-01 . 29188135E-02 . 15778093E-02 200 

6 0.47619048E-01 . 75109770E-02 . 28446236E-02 200 

7 0.47619048E-01 . 87010486E-01 . 36285359E-01 200 

8 0.47619048E-01 . 20008583E-01 . 10505322E-01 200 

9 0.47619048E-01 . 25432033E+00 . 12331680E+00 200 

10 0.47619048E-01 . 14796519E+00 . 51432260E-01 200 

11 0.47619048E-01 . 14443649E+00 . 47481858E-01 200 

12 0.47619048E-01 . 13431004E+00 . 57817377E-01 200 

13 0.47619048E-01 . 62963498E-01 . 18297859E-01 200 

14 0.47619048E-01 . 14057286E+00 . 83384476E-01 200 

15 0.47619048E-01 . 16790199E+02 . 67749 184E+00 200 

16 0.47619048E-01 . 61926860E-01 . 12581699E-01 200 

17 0.47619048E-01 . 12722393E+00 . 37044922E-01 200 

18 0.47619048E-01 . 15741819E+00 . 64646409E-01 200 

19 0.47619048E-01 . 16399987E+00 . 54156119E-01 200 

20 0.47619048E-01 . 21743058E+00 . 13237257E+00 200 

21 0.47619048E-01 . 13765610E+00 . 84254619E-01 200 



Born scanned with 4200 calls = . 19300621E+02 

The maximum value of the cross section= . 98788719E+03 

shifted to . 14818308E+04 



Weight 21 = 0.7132210911535855E-02 rescaled to . 7132210911536152E-02 
Calculating the Born cross section with the weights: 

1 0.12237203E-01 . 12237203E-01 

2 0.13250817E-01 . 25488019E-01 

3 0.46795867E-03 . 25955978E-01 

4 0.73460104E-02 . 33301989E-01 

5 0.15122900E-03 . 33453218E-01 

6 0.38915728E-03 . 33842375E-01 

7 0.45081704E-02 . 38350545E-01 

8 0.10366808E-02 . 39387226E-01 

9 0.13176795E-01 . 52564021E-01 

10 0.76663435E-02 . 60230365E-01 

11 0.74835153E-02 . 67713880E-01 

12 0.69588457E-02 . 74672726E-01 

13 0.32622525E-02 . 77934978E-01 

14 0.72833337E-02 . 85218312E-01 

15 0.86993051E+00 . 95514882E+00 

16 0.32085424E-02 . 95835737E+00 

17 0.65917015E-02 . 96494907E+00 

18 0.81561205E-02 . 97310519E+00 

19 0.84971292E-02 . 98160232E+00 

20 0.11265471E-01 . 99286779E+00 

21 0.71322109E-02 . 10000000E+01 



The maximum weight updated from . 14818308E+04 to . 16809386E+04 
The maximum weight updated from . 16809386E+04 to . 17190927E+04 



Final result : 

Integral = 0.206457E+02 Standard dev. 


= 0.18E+00 


with 


20000 


calls 


Efficiency of the events acceptance = 


0.129E-01 








Approximated cross section e+ e- -> t 


t -> 6f: 








Final result : 

Integral = 0.193058E+02 Standard dev. 


= 0.66E-01 


with 


20000 


calls 


Born cross sections (fb) : 


ecm cs sd cstta sda csbbww sdbbww 



0.5000E+03 0.20100E+02 0.27E+00 0.19306E+02 0.66E-01 0.20646E+02 0.18E+00 



